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Abstract 

Relativistic equation of state and velocity comparable with the 
speed of light are included in consideration of a superfluid rotating 
in a cylindrical container. Minimizing the free energy, we derive the 
equation of motion. It admits an analytic solution, the solid-body 
rotation inside and irrotational motion near the walls of the vessel, 
providing the vortex quantum is not extremely high that is satisfied 
for real astrophysical objects. The relativistic velocity of the vessel 
and the relativistic equation of state results to the deviation of the 
angular velocity of the solid-body motion inside and that of the vessel. 
The boundary between the solid-body and irrotational motion is also 
shifted sufficiently leading to a difference between the total angular 
momentum of superfluid and the normal matter. 



1 Introduction 

The purpose of the present study is to investigate the dynamics of the rela- 
tivistic rotating superfluid. It will be a generalization of the non-relativistic 
theory |1] . The medium is usually called relativistic in two senses: if it has a 
relativistic equation of state or when it flows at a relativistic velocity. Both 
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conditions (especially the first one) take place in neutron stars whose mas- 
sive interior is composed of a superfluid nuclear matter [[|. A superfluid 
with quantum vortices has a third measure of relativity; it corresponds to 
the vortices - whether they are relativistic ||. This possibility pertains bet- 
ter to the string fluids rather than to the real neutron star matter, although 
some small relativistic effects are also feasible here. We do not omit them 
immediately from the beginning and do not discuss more than it is necessary 
in the present paper with the aim to find possible applications to the neutron 
stars. 

The signature of the Minkowsky space with a metric tensor 
diag({ — h ++} and (if not specified otherwise) the natural system of units 
(H = c = 1) are used in the paper. 

2 A relativistic superfluid in a rotating cylin- 
der 

1. Let a superfluid is contained in a cylindrical vessel rotating along the axis 
x 3 with a constant angular velocity uo which is assumed to be large enough 
for existence of a vortex array in a superfluid, while the distance between 
neighboring vortices is small in comparison with the size of container Rq. 
The equilibrium condition is determined by a free energy minimization JTJ: 

5F = 5E-u i 5L i = Q (1) 

where 

E = J T ° d 3 R J = (0, 0, u) (2) 

is the energy, while 

Li = SijkLi L 3 = — L 12 (3) 

and H 

L jh = J (x j T k0 - x k T j0 ) d<f)RdR (4) 

is the angular momentum and its tensor form, respectively. The superfluid 
velocity V (R) at distance R from the axis of rotation follows immediately 
from the constraint 5F = 0. 
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r; = { p + p) u» Up - Wl w»°w ap - (7) 



2. In general the liquid and the vortices are not separated H due to the 
Lagrangian dependence A on the amplitude 

h 2 = h v h v (5) 

of the helicity vector 

h u = \e va ^^ a W M (6) 

which, as we see, is the cross product of chemical potential 1-form p a = 
pu a and vorticity 2-form Wp T If dA/dh = 0, the liquid and the vortices 
can be considered separately (as in the frames of the dilatonic model 0], 
which may be treated as an analogue of the two-constituent superffuid model 
U). Namely, the particle number current n v is collinear to the chemical 
potential flow: n u = nu u = $ 2 /i l/ . This situation occurs in the weak vorticity 
limit which takes place inside a typical neutron star ||. Thus the energy- 
momentum tensor of superffuid will be || 

A 

W 
where 

A = K& 2 $ 2 = - (8) 
A 4 

while pressure P of the superfluid constituent and its rest-mass density p 
satisfy the relation pri = p + P pertaining to a perfect fluid. The second 
term in the right side of (|7|) is negligible in the weak vorticity limit, while 
the pressure function \I/ — > P. 

3. For a vortex line situated at the distance R from the axis of the vessel 
it is convenient to switch to the reference frame rotating with the superfluid 
medium at the velocity V (R). The metric in the rotating frame has the form 

i 1 

ds 2 = - (l - oj 2 R 2 ) dt 2 - 2uRd<f dt + dR 2 + R 2 d V 2 + dz 2 (9) 

On account of small b (with respect to the size of the vessel) the metric in 
the vicinity O (R) of the vortex line can be regarded as flat. Thereby, we can 
also introduce the local cylindrical coordinates (t, r, <f), z) and calculate the 
invariant circulation integral |3j 

— f p v dx v = K = hn (10) 
2ir J 

dU 
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which determines the local superfluid velocity (r) round the vortex 

v f r ) = (±± = JL = 0k (ii) 

m mr r 

where m is the mass of the boson pair and l c = h/ (mc) is the Compton 
length. Note that the superfluid velocity (defined as a gradient of the bose- 
condensate wave function phase ||) does not coincide with the usual four- 
velocity 

\x v m 

u u = — = — v v (12) 
ix /i 

This solution (locally) satisfies the irrotationality condition (5|. || 

rot-u s = /io = const (13) 

(where 

ToW=Ld R (vR) (14) 
ti 

in cylindrical coordinates) which takes no place in the global sense, particu- 
larly on the axis of the vortex line. 

Although (in the local reference frame pinned to the vortex) the averaging 
over domain O (R) yields zero momentum 

M = o (15) 

the presence of a vortex results in the global rotation of the superfluid: 
V (R) 7^ 0. The world sheet £ of the vortex line plays the role of support of 
the vorticity 2-form W VQ (W ug ^ on E). Taking into account the link (|TT|) 
between the superfluid velocity v<p and the flow chemical potential /i^ and 
calculating the integral (ID) in the laboratory reference frame, we have 



\ r If 771 

k = — f fi u dx u = — / W ve do ve = —\rotv\ nb 2 (16) 

2,71 J Z71 J 2,71 

dU U 

where |rotu| (R) is the value averaged over domain U (R). Indeed, the global 
velocity field V (R) does not satisfy the irrotationality condition. Formula 
(|16D implies relation between the velocity field and the distance between the 
vortices: 

b2 = _2k 
m \Totv\ 
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Equations (jT0|)-(p~3l) also imply 



K 2 



p 2 = -fiffi e = £-'it (18) 



where /i* (R) = \//i /io = Po/y 9oo (R) and \i is invariant and can be treated 
as the rest-frame chemical potential. The four- velocity ([12|) constructed from 
(|18D will be 

tfr= J— , (-l,0,tP,0) «,( r ) = — = ^ (19) 
V 1 - w 2 \ir r 

where the modified Compton length is 

k , m , , 

r c = - = / c -n 20 
/i /i 

4. Followed the non-relativistic procedure ffl, we can take into account 
the possible non-uniform distribution of vortices if: calculate the energy of a 
single vortex line and multiply it on the density of the vortices 

N(R) = ± = ^^ (21) 
v ' nb 2 2ttk v ' 

Substituting (|19"D in (|7|) we find the energy 



E(R)= /T °rdr= / (f^ + fi-U (22) 



p (r) w 2 P 

1 — U> 2 1 — U? 2 

of a single vortex calculated in the local reference frame (associated with 
point R). Of course, due to the local velocity w (r) ^ round the vortex 
fll9|) the mass density p(r) deviates from the proper rest mass density p s . 
Taking into account the evident relation p{r) = p s \/l — w 2 between them 
and Eq. (^2|), we obtain 



E(R) 



71 



( / / \ 9 I b + \ b2 - r c\ 9 (b 2 -r 2 

Ps b^b 2 - r 2 - ay* a 2 - r 2 + p s r 2 In V + Pr 2 c In ' 



a + \ a 1 — ri 



a? — r 2 



(23) 

where a is the inner cutoff radius determined from the kinetics concepts. 
The typical values for neutron stars are M: b ~ 10~ 4 cm, while a is, at least, 



5 



several times greater than r c ~ 1 fm (we do not consider the exotic cluster 
vortices with n ~ 10 12 ). Thus, Eq. (p3|) is simplified so that 



£ (i?) = vr 



p s b 2 + p s r 2 c In 



2b 



+ Pr\ In 



a? — r 2 



a + y a 

Multiplying ( p4j) by (^1|), we find the superfluid energy density 



(24) 



e(i2) =p s + p s ^|ln 



at point R, or (more roughly) 
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a + \ a z — ri 



a 2 — r 2 C/ 



(25) 



e{R)= Ps + { Ps + 2P) r ^r- 



(26) 



It is the mean value (we may define it also as e s = (Tq)) of the energy density 
in the vicinity O (R), for (p5"D bears no dependence on the local coordinate r 
and, hence, does not reflect the fine structure of the vortex cell. The energy 
e measured in the laboratory reference frame will be given merely by formula 



1-V 2 



(27) 



where Eq. ( |T5D was taken into account and the global ordinary velocity Vr 
is resulted from the global superfluid velocity (do not mix it with the local 
filed ity) 

JS= (28) 

which coincides with Vr in the non-relativistic limit. 

The mass density p^ in the laboratory frame is expressed through mass 
density in the rotating frame as pn = PlV^- — V 2 , while pr is determined as 
the averaging over domain O (R): 



PR = (p(r)) 



1-V 2 



( r 2 b 
rdr ~ p s I 1 + In - 



(29) 
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Indeed, this term, which has appeared in includes contribution from 

the vortex. 
Thereby 

while, according to (|ll|), (0), and (|2 



where 



r 2 k m r c m y . _,. 
— = — rotf — = — rotf — = — rotf (31) 

9 

X = 2vr/m— (32) 

reduces to the usual circulation quantum |I[ in the non-relativistic limit. We 
shall bellow a more convenient expression 

X m nm 2 , N 

a = — = cr c — = cl c -— 33 
87T 4/x 4 fi 2 

instead. 

5. The total angular momentum density [ffl 

(0,0, L 3 ) ^ L 12 =M 12 + S 12 (34) 

includes spin contribution of the vortices S-* and the orbital part M^ k from 
the rotating fluid itself. The latter may be calculated (in the laboratory 
frame) by formula (|j): 



M 12 



/ x'T 20 dV = ( { PR + P) VrR d 3 R (35) 
J J Jl - V 2 



while the spin contribution from a single vortex is determined by integral 
(0): 

s 12 = — f fi u dx u = k (36) 

dO 

It is invariant with respect to transformation from the rotating to the lab- 
oratory reference frame. This is derived immediately from the constraints 
i if 

s% = s 2 p% = (37) 
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which are the same in any reference frame; here p M is the momentum. Since 
in the rotating frame 



-\ivp ' 



S "P 



(s°,0,0,s*) p»= (p°,0,/,0) 



(38) 



and p° 7^ 0, we find, by means of (|3"?D and metric (|9|), that s° = 0. Hence, 



(s 2 ) and, in the laboratory frame 



& P 



s°, 0, 0, s : 



= 0,^,0 



m 2 Miyp 
The identity (|37D then implies that 

because so = 0, according to ( flO|) , and the metric is defined as 
also s = and, hence, s 2 = s 2 = s 2 . 
Thus, we can define the spin density 

PR s n 
m 

measured in the rotating reference frame (related to the point R). 
Hence, 



(39) 

(40) 
Therefore, 

(41) 



S 
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i-vg 



<TR 



(42) 



will be the total spin in the laboratory frame || [7j. Combining Eqs. fl3~4D , 
fl35|), ( f42|) and (|29|), (|3~T|), we get the total angular momentum 



2np s 



RdR 

T^v 2 



b 2 a) p s 



V R R-[l + r hn-\^- 
\ br aim 



6. The total free energy is constructed from (]!]), (|30|), (|35D, 



(43) 
91) so 



2-nps 



Ps o 2 a 



b 2 a) p s 



ujR 



r 2 b\ kuj 
1 + -f In - — 

(44) 



RdR 

T^v 2 
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that is 
F 



J RdR Uj + ( 



„ X • -> * 7r|rott>|a 



2np s J { 
where, in the ligth of (|3l| ) 



vuR — S7) — \iotv\ In 

07T 



X 



7 



i-vg 



the parameters 



J 



m 



Ps 



m 



and the equation of state (EOS) index 



r = i + * 

Ps 

tend to a unit in the non-relativistic limit. Note that 

ku Tiuon L nuo c 



(45) 



m 



R2 



LU < Um 



Ri 



(46) 



(47) 



(48) 



(49) 



implies the ratio of the spin constituent of the rotation energy per particle 
(one may call it as "rotation quantum") to its rest-mass energy. It does not 
exceed the ratio of the Compton length to the radius of the star and occurs to 
be extremely small, namely, it is evaluated as 10 -21 even at angular velocity 
uj = 10 3 which is rather high for neutron stars. One may only mediate 
at application to a hydrodynamic model of rotating nuclei emphasized e.g. 
in [f|. 



3 The equation of motion 

Varying expression fl4"5|) over 8v, as shown in Appendix in details, we obtain 
the equation of motion 



J 



uRT+a Ey/l + v 2 - vuR 



dnWotv\ „ v 

■fCKC- 



rotf 



VT 



rotf 



v ea 

R lU ^ 

(50) 



—auR\rotv\ = 



9 



Or, including the speed of light explicitly 



dji\iotv\ 1 _ 

fa— p= I c^f — —In — 1 —at 




v r / l 

uoRT—J = = a E\l + v 2 /c 2 -vuR , , .,_ 

yjl + v 2 /c 2 L V c 2 J |ratu| c 2 yT+^/c 2 

(51) 

While its dimensionless form 

( J . = = -xv)w = ^ (g . V = (-\n^-w]w + xW 2 + (W- ZVTT^ 2 ) d x w) 
\ Vl + " / I VI + v 2 \x ¥ J v ' J 

(52) 

with 

i*=t"$ W = ±d x (xv)~\mtv\ x=f (53) 

is also convenient for further discussion. 

When simultaneously c — ► and r — > 1, Eq. (^) reduices to the well- 
known non-relativistic equation [|TJ 

v-uR + a d f Ot f =0 (54) 
| rotf | 

for the rotating superfluid helium. 

While the non-relativistic equation (|54| ) was describing the solid-body 
rotation v = ujR in the inner domain (R < Ri) and the irrotational motion 

roU7 = (55) 
in the outer domain (Ri < R < R 2 ), the general equation fl50D-([52"D also 



determines an irrotational solution (p5|) but it does admit any solid-body 
rotation in the strict sense. However, when the right side of (0) is small, the 
left side is splitted into a product of these two solutions which are connected 
in the transition region (R ~ R4) whose width I is small with respect to Ri 
and R2 — Ri- 

The intermediate region (between the solid-body rotation at small R and 
the irrotational motion in the outer layers) may become relatively broad at 
significant right side of (B2). On account of small £*, this may occur at ultra- 
relativistic velocities v ~ 1 /£ ^> 1 . Even if we apply the present analysis to 
the cluster vortices with n ~ 10 12 , the lowest value will be 

v > 10 7 (56) 
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Although one could consider this situation in the frames of the hydrodynamic 
nulear model, it does not appear when we study the real neutron stars. 
Therefore, except for the main transition region, the solution is determined 
by equation 

luRT 1 rotv = O (0 (57) 



VI + v 



2 



which yields 



V+ = TuR Ji + V l = (58) 
'l+v 2 + y/l - (TuR) 

or, the angular velocity 

n = Tlu (59) 

of the solid-body rotation in the inner region. It differs from the non- 
relativistic rotation by multiple T = 1 + y. So, for an ultrarelativistic 
equation of state (p = 3P) a superfluid rotates at angular velocity f2 = |u; 
rather than fi = u. For a stiff EOS (p = P) the superfluid angular velocity 
exceeds twice the angular velocity of the vessel. The conclusion applied to 
the neutron stars is evident: the angular velocity of the superfluid core Q 
differs sufficiently from the angular velocity of the crust u>. 
As for the irrotational solution 

v- = | (60) 



of (p?\) in the outer region, it must obey the boundary condition v (R2) = 

(61) 



uR 2 /Jl — oj 2 R\/c 2 = Q/R2 and, hence, be defined ultimately 

L0R2 R2 



That is 



- oo 2 R 2 2 /c 2 R 
xl q 



xJ 1 — Xq x 



(62) 



in the dimensionless form, where 

xl L0R2 

= x = — 

1-x 2 c 



q = r^— x = ^- (63) 



Indeed, in the light of formula (|58|) only xq < 1/T has physical sense. 
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4 The boundary between two solutions 



4.1 The extremum 

In order to find the boundary Ri between the two types of motion, deter- 
mined by formulae (^) and (|58|), respectively, we verify the free energy (|64l) 
over unknown Ri which must minimize F. The latter, on account of the 
thin intermediate region (whose size is determined by usual expression [|TJ 
multiplied by T), can be splitted into a sum 

R2 Ri 

F = J '/_ (R) RdR + J f+ (R) RdR (64) 

Ri R\ 

of irrotational /_ (R) = f [v- (R)] and solid-body contribution /+ (R) = 
f [v + (R)], which we present in convenient dimensionless form 



(65) 



+ 



2np s 




+ o(a)\ (66) 



where quantum contribution o (a) is, evidently, proportional to the quantum 
number n (the non-linear dependence appears at ultrarelativistic rotation, 
when the quantum term is sufficient, that may occur in nuclear hydrody- 
namics) . 

The condition of extremum 



dR 
or 



dF 

' R=Ri = -U(R t )R l + f_(Ri)Ri = (67) 



/+ ( Xi ) = /_ ( Xi ) (68) 

allows to determine radius Ri = Xic/u without direct calculation of the total 
free energy (|64"D, ([[]). 
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4.2 The non-relativistic velocities 



In principle, we may confine ourselves with the non-relativistic velocities of 
rotation (while the ratio Pj p is sufficient), because < ^ at typical values 
Q < 10 3 s _1 and R2 — 10 km for the majority of neutron stars. Even for 
very high fl, approaching to 10 4 s _1 [Kj, the quantity x 2 = | is small. 
At non-relativistic velocities the free energy (pi] ) reduces to 

F 

± nr 



, RdR TvcuR + Ealiotvlln— (69) 

2np s J 12 a 2 

and its irrotational and solid-body contribution are specified immediately: 

/_ (it!) = / [«_ (it!)] = _ r^i? (70) 

/ (it!) = / [u. (i2)l = + EaR\rotv\ In — - TQiuR 3 (71) 

2 1 a 2 



Substituting them in extremum condition (p7\) leads to 

2 R 2 OT1 2 fi 2 i? 2 2~a ool b 2 

^- 2r ^=-^ + ir 2ni v (72) 

and, after plain arithmetics, yields 



Vr V r 



The quantum term in the right side of (|73|) is at least several orders less 
than i? 2 . But for superfluid helium [|I|] (r — 1 = — ~ 10~ 16 ) it defines a 
narrow band 




i? 2 - = W - In — ~ 10 2 cm (74) 



of irrotational motion near the walls of the container. For a relativistic 
matter, whose pressure P is compared with its energy density p, the deviation 
of the radius 

«. S ^ (75) 

from i?2 is significant. For an ultrarelativistic matter Ri = 0.87i?2, and 
Ri = 0.71i?2 for a stiff matter (r = 2). Therefore, a 35 4- 65 % portion of the 
total volume (which is proportional to R 2 in cylindrical symmetry) is free of 
vortices. 
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4.3 The relativistic velocities 

Substituting fl6"5| ) and ( |66| ) in ( |68j ) we obtain the equation 



J i + 4 - rg = V 1 - r ^ 2 + («) ( 76 ) 

V -^j 

for dimensionless = uRi/c. As we have mentioned above, the quantum 
term o (a) does not play significant role in defining Ri at the relativistic 
equation of state. Solving ( |66|) without o(a), we find that 

i*. = q (f^£ - fj = (V4-4. ? + r^ - r*») (77) 

generalizes ( |75| ) at relativistic velocities. The ratio Xi/xo always decreases 
with the growth of Xo implying that relativistic rotation tends to diminish 
Ri. And the ratio 



•'■V 




(78) 



corresponds to the Yxq — > 1. Note, however, that Xi — > xo as soon as T — > 1; 
meanwhile, the smallest ratio X{/xq = 0.66 is achieved at T = 2. 



5 The angular momentum 

5.1 The general expression 

The total angular momentum (|43|) can be presented as a sum of classical and 
quantum terms 

— — = / RdRTvR + L(a) (79) 
2%p s J 

where the quantum contribution L (a) is negligible in comparison with the 
first term, because the right side of equation (|5^) is small (the quantum 
term is sufficient only at ultrarelativistic rotation). Therefore, omitting it and 
substituting the velocity fields flBTD and ( |58l) in (^), we get the dimensionless 
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expression 

XO 



L U _. /'_2 J _ Fa; , /"..2, 9 



x 2 cfa; = + / £ 2 cfe- (80) 



o 



2np s Tc 2 7 vi - r 2 x 2 / x 
easily integrated up to 



l ^ 2 - - r^ 2 (2 + r 2 * 2 ) i 2 

+ -g(a; -x i J (81) 



2?rp s Tc 2 3r 3 2 

where Xq = 00R2/C and determined by ( \T?() . 

5.2 The non-relativistic velocities 

The non-relativistic angular momentum can be derived either from for- 
mula (|5T) or from ( |7Sp Substituting there ( pT| ) and (0), we have 



2vrp 

that is 



J QRTujRRdR + J ^TuRRdR (82) 



Ri 



1 ,2 D 4 /- 1 



~-u 2 R\ T - - (83) 
2vrp 2 2 V 2/ V ; 

and always exceeds the non-relativistic value ^ = \u 2 R\. Their ratio is 
5/3 for T = 4/3 and 3 for T = 2. Formula (j83|) may be also compared with 
the angular momentum 

^- = \u 2 B*T s (84) 
2vrp 4 2 V ' 

of normal fluid rotating as a solid body with the same radius R2 and angular 
velocity u. If normal matter rotating with angular velocity u n transfers to 
superfluid with the same angular momentum, its angular velocity 



u 2 T 2 



loI 2Y-1 

will be always greater than u n . The first relativistic correction 

u 2 s r 2 r 3 f„ at - 3 \ x 2 



(85) 



+ -!— r- a + oM) (86) 

2r - 1 1 (2r- 1) 3/2 / 3 v °y v ; 



c^ 2 2r - 1 2r - 1 V ( 2 r 
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does not contradict to ( p5|) states the same, namely that the superfluid ro- 
tates also faster than the normal fluid with the same angular momentum. 
Meanwile, 



<4 



(87) 



when r — > 1. Indeed, the difference between oj s and u n is not negligible at 
relativistic EOS, because T — 1 is sufficient. This forms the background for 
further application to neutron stars and pulsar glitches. 

5.3 The relativistic velocities 

The angular momentum of relativistic superfluid (|8~l|) differs from its expres- 
sion 



L n _ 2 - y/1 - (2 + T*xp 1 

2^rv " 3f3 x ° - f (88) 

for normal matter with the same parameters, i.e. density p s , equation of 
state (indicates T), external radius (dimensionless xq). The angular momen- 
tum (|88|) , as well as (|8l|) , inreases with growth of Xg. While always L s < L n , 
especially at relativistic rotation. Indeed, they coincide at T — 1 (for we 
have omitted the negligible quntum contribution), while their ratio (L s /L n ) 
at different T is given in the table below. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Tx 
L s /L n 

T = 1.1 0.992 0.99 0.989 0.987 0.984 0.979 0.973 0.964 0.947 0.911 0.694 

T = 1.2 0.972 0.968 0.964 0.958 0.95 0.94 0.926 0.907 0.876 0.82 0.60 

T = 1.6 0.86 0.85 0.83 0.82 0.80 0.78 0.75 0.71 0.67 0.60 0.42 

T = 2 0.75 0.735 0.72 0.70 0.68 0.65 0.62 0.59 0.54 0.48 0.33 



Instead of x , the universal variable Tx stands here in the first row. 



(89) 



6 Conclusion 

Summarizing, we list the main results. Having derived the equation of mo- 
tion (|5Tf )-(|52"l) for the rotating relativistic superfluid, we find that it may be 
reduced to a product of two independent solutions (|5T|), namely the solid body 
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rotation fl58| ) in the inner domain R < R4 and the irrotational motion (|6T[ ) in 
the outer domain R > R4 with a relatively thin region where they interfere 
at R = Ri. The splitting (0) into these independent solutions was possible 
due to a small rotation quantum (|49|) implying that the velocity of rotation 
is not extremely high (|56D , that takes place for all real macroscopic objects, 
particularly, the neutron stars. Although, considering the nuclei within the 
hydrodynamic approach, one encounters the opposite situation and has to 
solve equation of motion (|52|) numerically. 

As for solutions (|58D and(^l|), there should be noted the fact of depen- 
dence on the equation of state, besides the relativistic expression for velocity. 
The relativistic equation of state and relativistic rotation are characterized by 
dimensionless quantities T = 1 + P/ p s and x$ = ujR^/c, respectively (which, 
in the non-relativistic limit, become r — ^ 1 and Xq — ■> 0). However, without 
regard of rotation rate the relativistic EOS results to difference between the 
angular velocity fl5"9|) of the superffuid - obeying the solid-body rotation (|5g|) 
within the vessel, - and the angular velocity uj of the vessel itself. In other 
words, the superfluid rotating as a solid body rotates at angular velocity 
which T times higher than that of the vessel. 

The same quantity V appears in formula (|73|), (|76|) for radius Ri which 
corresponds to the boundary between the solid-body and irrotational mo- 
tion. While the non-relativistic value of Ri deviates slightly from R2 due to 
the presence of quantum term (|74|), this quantum term plays an insufficient 
role when the difference (|75"|), (\n\) between Ri and R2 becomes consider- 
able on account of the relativistic equation of state or relativistic rotation. 
Particularly, the boundary radius Ri is determined merely by T for a rela- 
tivistic matter rotating at low velocity ([751). Or, briefly, Eq. (|77|) states that 
both relativistic rotation and EOS shorten the distance Ri, i.e. the region of 
irrotation motion is much wider than that in non-relativistic helium. 

So, the quantum term involved in the relevant non-relativistic formulae 
becomes insufficient for a relativistic superfluid. Thus, the total angular 
momentum (|79|) is given approximately by formula ([81]). The angular mo- 
mentum of a slowly rotating relativistic superfluid ( |8BD is always less than the 
angular momentum ([84]) of the normal fluid rotating with the same angular 
frequency. Hence, if the normal matter, rotating with frequency u n , trans- 
fers into a superfluid state (with angular momentum conserved), its angular 
velocity u> s will be increased (^5[). The general dependence (|8l| ) of oJ s /u n 
on the initial angular velocity xq = 00R2/C and the EOS index T is more 
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complicated as illustrated in table (B^). However, the angular momentum of 



superfluid (81) occurs to be always smaller than that of normal matter (j38|) 
and the relativistic rotation contributes to this tendency. 

Thus, the equation of state and velocity uR 2 of a rotating superfluid are 
of particular importance when they belong to a relativistic range. A plenty 
of applications (the pulsar glitches, for instance) are expected to be derived 
from here. And further development in the light of nuclear hydrodynamics or 
solution with a solid core (that is R2 > R > R\ > 0) may be also proposed. 



7 Appendix: variation of the free energy 

Since the rotation quantum (03) is small for usual macroscopic objects, we 



can rewrite (f45|) as 



F 

Ps 



2 

2tt J f RdR f = 7 J + a|roti7| In ^ {vtoR - ^7} - TvuR (90) 



Firstly we specify variation 



BC 

SG [v] = —Sv = G' Sv 
ov 



(91) 



of an arbitrary function G (v). Hence, the simplest expressions for variations 

' Jv 



5B = 5(-fJ- TvuR) 



for 



7 



i-vg 



vTT 



TtuR Sv 



The variation of 

A = —a\rotv\ G In 
is performed so: 
5 A = a 



X d R {vR) a 2 

— 1 ^7 = a 7^ G m — 

7T roto \a l R tr 



S(l/b 2 



d R (5vR)G\n- + d R (vR) ]/j)2 



■G + d R (vR) In — 8G 



(92) 
(93) 

(94) 



dR = 5A 1 +5A 2 +SA 3 
(95) 
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where 



In 



b 2 



1/b 2 



d R \i:otv\ 


5 (1/b 2 ) 


S\iotv\ 


d R (SvR) 


\rotv\ 


1/b 2 


\rotv\ 


d R (vR) 



(96) 



Hence, 



(97) 



2 

SA 2 = J d R (SvR) GdR 5A 3 = J d R (vR) In ^G' 5vdR 
The first two terms SA\ + 5A 2 of (pq ) are simplified as 

2 2 / 2 \ 

<L4i+<L4 2 = y d R (SvR) In ^-Gdfl = SvRln ^G\f-j 8vRd R f In ^-Gj di? 

(98) 



The total variation of ( |95| ) will be 

,2 



<5A = tfwiZln ^-G|g + J Sv In ^ G' [v] - itf* (g [v] In ^-j | dR 

(99) 

Taking into account that variation 5v vanishes at the edges of the vessel (i.e. 
at R = Ri and R = R 2 ), and adding (j92|) to the latter formula, we find the 
total variation of free energy 



8F 



Sv 



Jv 



VTT 



TuR + G' (v In i- - Rd R v) - Gfl ^ffl l rfi? 
\ cr y |rotv| 



(100) 

which for arbitrary <5t> and G [v] = vujR — ipl [v] yields equation of mo- 
tion (0). The relevant non-relativistic version is immediately obtained if we 
put G = 1. 
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